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Analytical multi-domain solutions to the dynamical (Landau-Lifshitz-Gilbert) equation of a one- 
dimensional ferromagnet including an external magnetic field and spin-polarized electric current are 
found using the Hirota bilinearization method. A standard approach to solve the Landau-Lifshitz 
equation (without the Gilbert term) is modified in order to treat the dissipative dynamics. I estab- 
lish the relations between the spin interaction parameters (the constants of exchange, anisotropy, 
dissipation, external-field intensity, and electric-current intensity) and the domain- wall parameters 
(width and velocity) and compare them to the results of the Walker approximation and micromag- 
netic simulations. The domain-wall motion driven by a longitudinal external field is analyzed with 
especial relevance to the field-induced collision of two domain walls. I determine the result of such 
a collision (which is found to be the elastic one) on the domain-wall parameters below and above 
the Walker breakdown (in weak- and strong-field regimes). Single-domain- wall dynamics in the 
presence of an external transverse field is studied with relevance to the challenge of increasing the 
domain-wall velocity below the breakdown. 
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I. INTRODUCTION 

Description of the magnetic-field- and electric-current- 
induced motions of domain walls (DWs) in nanowires has 
became a hot topic because of novel methods of storing 
and switching the (magnetically encoded) binary infor- 
mation. These proposals offer a progress in the minia- 
turization of memory and logic elements, utilizing crucial 
advantages of magnetic information encoding, (when a 
bit is identified with a single magnetic domain) . Such in- 
formation is insensitive to the voltage fluctuations while 
its maintenance does not cost any energy, which enables 
the data processing with the production of a small heat 
amount. Currently investigated random-access memories 
are built of metallic nanowires, formed into a parallel- 
column structure, which store magnetic domains sepa- 
rated by DWs. Such a three-dimensional (3D) magnetic 
system has the potential of storing more information than 
devices based on 2D systems, like hard-disk drives or 
electronic memories, in a given volume 0,1. Also, an 
interesting concept of logical operation via transmitting 
magnetic DWs through nanowires of specific geometries 
is being developed [3|, I mention that ferroelectric 
nanosystems offer similar capabilities while their basic 
properties are studied with the same dynamical (Landau- 
Lifshitz-Gilbert) equation even though the effects of elec- 
troelastic coupling are strong [j| @ . 

In order to write and switch information, one can move 
the DWs via the application of an external magnetic 
field (parallel to the easy-axis) or via the application of 
a voltage which induces the spin-polarized electric cur- 
rent through the DW. The directions of the field-driven 
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motion arc different for the tail-to-tail and head-to-head 
DWs; thus the magnetic field induces the DW collisions, 
while the direction of the voltage-driven motion uniquely 
corresponds to the current direction. Field-driven mo- 
tion and current-driven (below the Walker breakdown) 
motion are possible due to the magnetic dissipation and 
its description demands inclusion of the Gilbert term into 
the Landau-Lifshitz (LL) equation. However, existing 
many-domain analytical solutions to the LL equation do 
not include the dissipation jj, 8], while approximate so- 
lutions using the Walker ansatz describe a single DW 
only Q. Since the parameters of the DW solutions to 
the Landau-Lifshitz-Gilbert (LLG) equation determine 
accessible values of technological characteristics (e.g. the 
minimal domain length, the bit-switching time, etc.), it 
is of interest to know the analytic DW solution to the 
LLG equation. Knowledge of the many-domain solution 
is of importance for preventing unwanted DW collisions 
which can result in an instability of the record and it en- 
ables verification of DW-collision simulations regardless 
of the internal structure or the geometry of the simulated 
system. 

In the present paper, I perform an analytical study of 
the dynamics of multi-domain systems including the dis- 
sipation. The dynamical LLG system is bilinearized fol- 
lowing the Hirota method of solving nonlinear equations 
and it is extended, via doubling the number of freedom 
degrees and the number of equations, into a time-reversal 
invariant form. The field solving the extended system 
contains proper and virtual (unphysical) dynamical vari- 
ables and the physical components of the solution are 
shown to satisfy the pimary LLG system in a relevant 
time regime. In particular, aiming to analytically de- 
scribe the field-induced collision, I establish asymptotic 
three-domain magnetization profiles (relevant in the time 
limits t — ¥ ±oo). With connection to the phenomenon of 
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the Walker breakdown, (a cusp in the dependence of the 
DW velocity on the external field and in its dependence 
on the current intensity), I modify the LLG model in a 
way to make it applicable below the breakdown. I re- 
duce it into a model of plane rotators. In this weak-field 
regime, the spin alignment in the DW area is saturated 
to a plane while, above the breakdown (the strong-field 
regime), the spins rotate about the magnetic-field axis 
(the easy axis). These considerations supplement micro- 
magnetic simulations of the DW collisions in terms of the 
studied DW-parameter regimes [l(|. Within the present 
method, I verify the Walker-ansatz predictions on the 
current-driven DW motion above and below the break- 
down, (including adiabatic and non-adiabatic parts of the 
spin-transfer torque) jllj|-fl5j. thus, showing the applica- 
bility of the present formalism to the field-driven motion 
of multi-domain systems with an additional voltage ap- 
plied. With relevance to the challenge of controlling the 
maximum DW velocity below the Walker breakdown, we 
analyze the longitudinal-field-driven motion and current- 
driven motion of the DWs in the presence of an additional 
perpendicular (with respect to the easy axis) field. 

A complementary study of the DW collision of mag- 
netic DWs in the subcritical regime is performed in a 
separate paper (l6j . 

In Sec. II, I extend the dissipative equations of motion 
of the ferromagnet in such a way as to make equations 
applicable to the unlimited range of time t 6 (—00,00). 
Section III is devoted to the analysis of its single- and 
double-DW solutions in the presence of a magnetic field 
and an electric current. I study the field-induced collision 
in detail. The plane-rotator approach to the DW dynam- 
ics is described in Sec. IV. In Sec. V, consequences of 
the application of an external field perpendicular to the 
easy axis for the DW statics and dynamics are consid- 
ered. Conclusions are given in Sec. VI. 



II. DYNAMICAL EQUATIONS 

The dynamics of the magnetization vector m (|m| = 
M) in the ID ferromagnet is described with the LLG 
equation 

d 2 m j3i , ~. ~ 

— + 7m x H + — (m • i)m x 1 
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The first term of the right-hand side (RHS) of (fJJ re- 
lates to the exchange spin interaction while the second 
term depends on the external magnetic field H; thus, 7 
denotes the giromagnetic factor (up to its sign). The 
constant /?i(2) determines the strength of the easy axis 
(plane) anisotropy and i = (1,0,0), j = (0,1,0). Note 
that the long axis of the magnetic nanowire is an easy 
axis for the majority of real systems; however, another 



choice of the anisotropy axes does not influence the mag- 
netization dynamics. The constant 5 is proportional to 
the intensity of the electric current through the wire and 
6 changes its sign under time-arrow reversal (the inver- 
sion of the electron flow) (T3, HH ■ The non-adiabatic part 
of the current-induced torque (which depends on f3) is of 
dissipative origin; thus, (3 — > with decreasing Gilbert 
damping constant a — > (one takes a, f3 <C 1). Its inclu- 
sion is necessary if one describes an observed monotonous 
motion of the DW below the Walker breakdown [Tl"j]-[l4j]. 
Notice that including the magnetic-dissipation term fol- 
lowing the original LL approach (changing the last term 
of JI]) into —am x [m x h c g], where h e &j — —STi/Srrij, 
and TL denotes the Hamiltonian) would lead to changing 
the constant (3 into (3 — a, [13J. Although the discussion 
of the relevance of both approaches to the magnetic dis- 
sipation remains open fljj , fl9| . I believe, the clinching 
argument for the Gilbert approach is the expectation for 
the proper dissipative term to be dependent on the time 
derivative of the dynamical parameter m. In other case, 
the dissipative term could influence static solutions to 
the LL equation while one expects the magnetic friction 
to be the kinetic. 

Since (QJ is valid only when the constraint |m| = M is 
satisfied, I intend to write equations of the unconstrained 
dynamics equivalent to ([T]). Introducing the complex dy- 
namical parameters m± — m y ± im z , I represent the 
magnetization components using a pair of complex func- 
tions g(x,t), f(x,t). This way I reduce the number of 
independent degrees of freedom. The relation between 
the primary and secondary dynamical variables 
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/Vfl + s*// f*/g + g*/f 

[where (•)* denotes the complex conjugate (c.c.)] ensures 
that |m| = M while there are no constraints on g, f. The 
transform enables bilinearization ("trilinearization") 
of (TT]) following the Hirota method of solving nonlinear 
equations [7], |8j- In the particular case H y = H z = 0, 
from (TT|) and @, we arrive at the trilinear equations for 

/, 9 

f [-iD t + JD 2 X + S(f3 - i)D x + aD t ] f* ■ 9 



+Jg*D 2 x g .g-\ 1 H x + p l 



\fV9 



g* [-iD t - JD 2 X + 5(p - i)D x + aD t ] f* ■ 9 

-JfD 2 J* ■ f* + (~7H X + A + y) isi 2 /* 

+yV/ = 0, (3) 

where D t , D x denote Hirota operators of differentiation 
which are defined by 

D?DZb{x,t) ■ c(x,t) = (d/dt - d/dt) m 
x(d/dx-d/dx'rb(x,t)c(x',t')\ x=x , t=t ,. 
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The inclusion of the dissipation into the LLG equation 
is connected to breaking the symmetry with relevance 
to the time reversal. Therefore, neither ([TJ nor ([3]) can 
describe the magnetization evolution on the whole time 
axis. In particular, the application of an external mag- 
netic field or current to the DW system (or creation of 
a domain in the presence of an external field) initiates a 
nonequilibrium process of the DW motion. Such a mo- 
tion cannot be present in the distant past (t — > — oo) since 
a nonzero value of the dissipative function [relevant to the 
Gilbert term in ([T])] would indicate unlimited growth of 
the energy with t — > — oo. Thus, solitary-wave solutions 
to (QJ) are relevant only in the limit of large positive val- 
ues of time fJl . Udi L This fact makes impossible an exact 
analysis of the DW collisions using ([3]) and motivates 
extension of the dynamical system within a formalism 
applicable to the whole length of the time axis. 

We write modified equations of motion using a similar 
trick to the one proposed by Bateman with application 
to the Lagrangian description of the damped harmonic 
oscillator [20(. It is connected to the concept by Laksh- 



manan and Nakamura of removing the dissipative term 
from the evolution equation of ferromagnets via multiply- 
ing the time variable by a complex constant pTj . how- 
ever, it demands an improvement in the spirit of Bate- 
man's idea [221 ] . The concept is to extend the dynami- 
cal system doubling the number of freedom degrees and 
adding the equations which differ from the original ones 
by the sign of the dissipation constants. Resulting ex- 
tended system is symmetric with relevance to the time- 
arrow reversal, however, its solution consists of physical 
and virtual fields. Let us mention that different quantum 
dissipative formalisms (non-equilibrium Green functions, 
thermo-field dynamics, rigged Hilbert space) are based 
on the Bateman's trick [231 ]. 

We extend the system of secondary dynamical equa- 
tions ([3]) since it describes unconstrained dynamics un- 
like the primary LLG equation JT]). We replace g, g* , /, 
/* in ^ with novel fields of the corresponding set g\ , g 2 , 
f2, fi and of the set of their c.c. For a,(3 — 0, in the 
absence of dissipation, 31 = 32 = 3, A = A = A For the 
case H y = H z = 0, the secondary dynamical equations 
transform into 



f 2 [-iD t + JD 2 X + 5((3 - i)D x + aD t ] f* ■ 9l + Jg* 2 D 2 x9l ■ 9l - (^H x + ft + y) /a/*0i - y/J/Tff* = 0, 
g* 2 [-\D t - JDl + S((3 - \)D X + aD t ] f* • g x - Jf 2 D 2 x f* • /* + (- 7 H X + ft + y ) 9* 2 9ifi + y<?a<7iA = 0. (4) 

Writing ((4]), we have replaced the last terms on the lhs of ([3]) in a way to be linear in g% > /a* > which ensures that 
they vanish (diverge) with time in the presence of H x 7^ with similar damping (exploding) rates as all other terms 
of these equations, (in particular, their damping does not modify the anisotropy). The additional equations of the 
dynamical system differ from Q by the sign of the dissipation constants a, j3 

h [-iD t + JDl + 5{-P - i)D x - aD t ] /I • 92 + JglD 2 x g 2 ■ g 2 - f^H x + ft + y ) /i/|<?2 - yAVafla = 0, 

g{ [-W t - JDl + 6(~(3 - i)D x - aD t ] f 2 ■ g 2 - J.hD 2 x f 2 ■ f 2 + f- 7 H x + ft + glg 2 fZ + ^3132/2 = 0. (5) 



Though the previous dynamical variables g{f), g*{f*) 
were mutually independent, they had to be c.c. to each 
other in order that the system of equations Q and their 
c.c. was closed. In the system of eight equations; ((4])- 
([5]) and their c.c, gx(fi) is not a c.c. to 32(72), while 
comparing (|4]) and ((5|), one sees that g 2 (x,t) (f 2 (x,t)) 
can be obtained from gi{x,t) {f\{x,t)) via changing the 
sign of its parameters a, ft 

Under the time-arrow inversion, the system of the 
novel equations transforms into itself if one accompanies 
this operation by the transform of the novel dynamical 
variables g 1[2) — > / 2 (i), A(2) -> -32(1) • The equations 
(|4]) and their c.c, which determine the magnetization 
dynamics for large positive values of time (in particular, 
for t — > 00), contain the differentials of the functions gi, 
9ii /l) fi ■ Therefore, the magnetization vector should 



be expressed with these functions in the relevant time 
regime. Writing the magnetization in the form 
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- (6) 



ensures that their components satisfy |m| = M, m x = 
m*, and they reproduce © for a — /3 = 0. In the regime 
of large negative values of time, in particular, for t — > 
—00, we can analyze the evolution of the magnetization 
with the inversed time arrow. It is described with the 
reversed magnetization vector 
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III. DOMAIN-WALL MOTION 

Let us analyze the multi-domain solutions to ©-© in 
the absence of any external magnetic field, H = 0. We 
search for the solutions which describe a single DW and 
two DWs in the forms /* = 1, g\ = w\e kl ' x ~ llt , and /* = 



1 + v * e kix-ht e k 2 x-i 2 t^ gi = WiQ . 

spectively, where kj = Re(fc J ), and sign(fci) = 
Inserting these ansatz into ((I])-©, one finds 



-^kix — lit 



w 2 e k2X - ht , re- 



-sign(fc 2 ). 
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j-^- (Jk 2 -fa) [Jk] - (fa + fa)] 



+5{l + i(3)k ] \, 



kj — 



lfa+fa(w 3 +w*) 2 /(4\ W j\ 2 ) 



J 



(8) 



The single- wall (two-domain) solutions with \k\\ 6 
(y/fajj, y/{fa + fa ) / J) describe moving solitary waves 
(topological solitons), @, H3- One sees the correspon- 
dence between the wall width and the spin deviation 
from the easy plain since the dependence of kj on Wj. 
When 5 — 0, static t wo-dom ain solutions represent the 
Bloch DW, (| fci| = y/faJJ, Wx = -w\), or Neel DW, 
(| fci| = yf(fa + fa)/ J, w\ — w*), respectively (the 
nomenclature of [24[ which differs from Neel and Bloch 
wall definition of e.g. [26|). These two solutions corre- 
spond to the ones found in (2f| [27j within the XY model. 
The DW profiles can be described with the functions 

-ilmiif 



rn+(x, t) = M — — — : sechffcix — ReZii + loglwi I 

HI 

m x (x,t) = — JVftanh[fcia; — Reiit + log|iui|], 



(9) 



(ki > relates to the head-to-head structure while k\ < 
to the tail-to-tail structure). At the time points of the 
discrete set t = nn/Ivah, n = 0, ±1,±2, one finds 
.9i = 92 = 9, fi = h = / an d © coincide with ©. 
Therefore, ([9]) is a solitary- wave solution to ([T]), (in par- 
ticular, it coincides with the one representing a Bloch 
DW or a Neel Dw for 5^0, [12]). Throughout the paper, 
we focus our attention on the externally-driven dynam- 
ics of the Bloch and Neel DWs, since these initially-static 
structures are the most important with relevance to the 
magnetic data storage. 

We establish that static double- wall (three-domain) so- 
lutions to the LLG equation cannot be written with the 
above Hirota expansion when k\ = — k 2 . In this case the 
coefficient v, 



fajk 2 wlw 2 



(Jk 2 - fa)(Jk 2 - fa - fa) 



(10) 



diverges with |fci| -> y/fa/J or \ki\ -> yj {fa + fa)/ J- 
Analogously to the XY model, the Hirota expansion 
is inapplicable to static three-domain configurations of 
the Bloch walls or Neel walls, while there exists static 



solution to ©-© which describes a pair of different- 
type (Neel and Bloch) walls (28|. In particular, for 

h = y/fa/J, k 2 = ~y/(fa + fa) / J, and W\ = 
w 2 = u>| , one finds 

faw^ 



2fa+fa~2y/fa(fa+fa) 



(11) 



Let us emphasize that we have not excluded the coex- 
istence of a pair of Neel walls or Bloch walls in a mag- 
netic wire. However, the overlap of both the topological 
solitons induces their interaction which leads to an insta- 
bility of their parameters and, unlike for nontopological 
solitons, is not a temporal one (29l |. 

Solving ©-© in the presence of a longitudinal mag- 
netic field H x ^ 0, we apply the ansatz 

f* = (l + v * e k ^-ht e k 2 x-l 2 t^ e ~fH x t/ (-2i+2a) ^ 

,9i = ( Wl e klX - llt + w 2 e k2X ~ l2t ) e -^ t /("2i+2a) ( 12 ) 

at the discrete time points t — t„ = 47rn(l + a 2 )/(^H x ), 
where n — 0, ±1, ±2 . . ., (let 5 — for simplicity). Behind 
these time points, in the presence of the longitudinal field, 
the last terms on the lhs of ©-© change faster (they 
oscillate with the three-times higher frequency) than the 
other ones. Therefore, taking the above ansatz, we apply 
an approach similar to the 'rotating wave approxima- 
tion' in the quantum optics. Since this ansatz describes 
the spin structure rotation about the x-axis, it is appli- 
cable when the external field exceeds Walker-breakdown 
critical value, \H X \ > Hw- From the single- wall solu- 
tion, (the case of w 2 — or w\ =0), for fci = y/fa/J, 
k 2 = —yj{fa + fa)/ J, we establish that applying the 
magnetic field in the easy-axis direction drives the DW 
motion with the velocity 



ci(2) = l\H x \a/[\k 1{2) \(l + a 2 )] 



(13) 



Correspondingly, for H x = 0, applying the electric cur- 
rent through the initially static wall drives it to move 
with the velocity 



5(1 + afa 



1 



(14) 



which is independent of the DW width. 

The essential difference between both kinds of the 
driven motion emerges from the analysis of the three- 
domain solutions. Under the external field, the two con- 
secutive DWs move in the opposite directions. The walls 
which are closing up to each other collide and eventually 
they can annihilate or wander off each other. The ap- 
plication of the electric current along the magnetic wire 
drives both the DWs to move in the same direction with 
the same velocity. Analyzing long-time limits of the mag- 
netization vector in different regions of the coordinate x, 
we establish the consequences of the field-induced colli- 
sion of the complex of a Bloch DW interacting with a 
Neel DW. We use the ansatz (TT21 and assume 5 = 0. 

Let i]j = kj(x — x j) — ~/H x at/(l + a 2 ), fjj = kj(x — 
x j) + jH x at/(l + a 2 ). For H x > 0, at t — t n (within 
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the above 'rotating wave approximation'), we find the 
distant-future limit of the magnetization (JB]) 
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-ijH x t/(l+a 2 ) 



1 + M 2 / 
1 



\wj\ 2 e 2f i" 



-M 



1 + HVIwjpe 2 ^ ' 



(16) 



and j ^ k. Identifying the parameters xoj with the 
DW-center positions, we introduce the restriction on Wj, 

\v\/\ Wj \ = l. 



,(2) 



We notice that rn^P , m^p as well as ni^\ 



are the Walker single-DW solutions to the primary 
LLG equation which describe the motion of well sepa- 
rated DWs [ij |3(J. Thus, our three-domain profiles of 
the fields ^ tend to satisfy (JTJ) in the limit t — » 00 ac- 
cording to the requirement formulated in the previous 
section. 

In the distant-past limit, we describe the magnetiza- 
tion evolution with the reversed time arrow. Following 
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where 
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u y , M ,* e »)3 e -i7^t/( 1 + a2 ) 

i + M7M 2 e 2l b- ' 



ml 3) = M 



1 



2 /\w k 
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1 + M 2 /M 2 e 2 * 



(18) 



and j ^ k. In order to consider the collision of the pair of 
DWs which are infinitely distant from each other at the 
beginning of their evolution, we determine the magneti- 
zation dynamics in the limit t — > — 00. For this aim, one 
has to invert the propagation direction of the kinks of rh 
and to reverse the arrow's head of the field vector m. Uti- 



lizing the properties to, 
(x + x ok ,0) 



TOa; 



mi j) 



(-x + 



£ fc,0)=TO + (-X + XQk, 0), 

^ofejO), we arrive at 



m+ix, t) 
m x (x,t) 



-ifvt\-x + 2xoi »*) Vi > ? /2 ~ 

-m\ (-x + 2xo2,t) 772 > m. ~ 

-TOa'^-X + 2x l,t) ?7i > T] 2 ~ 

-m^ (-x + 2.X02, t) 772 > 771 ~ 



(19) 



The applicability of the above procedure to study the 
asymptotic evolution of a single DW is easy to verify 
since any single-DW solution satisfies 



(20) 



m+(x, t) = — to+(— x + 2xoi,7j), 
m x (x, t) = fh x (—x + 2x Q i,t). 



Typically, one should consider the formulas ([15)) . ([IT 
with relevance to the case /?i ^> P2, thus k\ ~ — fc 2 , which 
corresponds to commonly studied crystalline magnetic 
nanowires, e.g. for Fe, FePt, /3 2 //3i ~ 10 _1 , [3l|. For 
noncrystalline (permalloy) nanowires Fei-^Ni^ deposited 
on a crystalline substrate, the easy-axis anisotropy con- 
stant determined from uniform-resonance measurements 
was found to be, unexpectedly, as big as in the crystalline 
nanowires [32j. Therefore, even when neglect structural 
effects in real systems which lead to the saturation of 
the spin alignment in the DW area to the easy plain 
or hard plain (the Walker breakdown) which suppresses 
their spontaneous motion, the spectrum of spontaneously 
propagating DWs in nanowires would be very narrow and 
their velocities would be very small. 

According to (fTS"]) , (JT5J) , two initially closing up DWs 
have to diverge after the collision. If one of the colliding 
DWs that was initially, for t — » — 00, described with the 

it is finally, for t — > 00, de- 



field ingredient fr$) , 



~ C?) 

TO 2; , 



U) 



Therefore, 



scribed with the field ingredient to 
reflecting DWs exchange their parameters xqi <-> £02, 
wi o w 2 , fci o — fc 2 . It is connected to exchanging the 
directions of the spin orientation in the yz-plane in the 
wall areas, (after the collision the Neel wall changes into 
the Bloch wall and vice versa as shown in Fig. 1). Our 
prediction corresponds to the result of the collision anal- 
ysis performed for spontaneously propagating DWs (in 
absence of external field, electric current, and dissipa- 
tion). According to findings of [13, HI], the DWs reflect 
during the collision in a way that one can say they pass 
through each other without changing their widths and 
velocities, however, with changing their character from 
the head-to-head one into the tail-to-tail one and vice 
versa. 

Up till now, we have considered systems of infinite do- 
mains whose energy cannot be defined. However, the 
smaller a domain is the bigger percentage of the Zeeman 
part of its energy is lost per time unit due to the DW mo- 
tion. The condition of the domain-energy minimization 
determines the direction of this motion. The domains 
aligned parallelly to the external field grow while the do- 
mains aligned antiparallelly to the field diminish. Any 
DW reflection induces a motion which contradicts this 
rule. Such a motion has to be decelerated and, eventu- 
ally, it has to be suppressed when the decrease of the DW 
interaction energy equals the increase of the Zeeman en- 
ergy. The outcome of a many-collision process in a finite- 
size system is the appearance of a ID magnetic-bubble 
structure similar to widely known 2D bubble structures 
[34l ]. The bubble size and concentration depend on the 
magnetic field intensity. Each bubble is ended with a 
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FIG. 1. The magnetization dynamics of a system of one 
Neel DW and one Bloch DW in a longitudinal field above 
the Walker breakdown. Their reflection takes place in the 
time region (At, 2 At) and it is accompanied by a change of 
the Bloch wall (of the velocity 112) into the Neel wall (of the 
velocity vi) and vice versa. Since \H X \ > Hw, the spin struc- 
ture monotonously rotates about the a;-axis. 



Neel DW at one of its sides and with a Bloch DW at the 
other side. Some analogy with a complex boundary of 
hard (quasi-2D) magnetic bubbles can be noticed since 
such a border contains alternating Neel and Bloch points 
in its structure [35] . Let us mention, that interesting con- 
cepts of storing and transforming the binary information 
have had been developed several decades ago with rele- 
vance to 2D magnetic bubble systems, though, they have 
been abandoned because of technological problems of the 
time (36j . 

Numerical analyses of the field-induced DW collision 
(micromagnetic studies) have been performed with rele- 
vance to flattened nanowires (quasi ID nanostripes) be- 
low and just above the Walker breakdown using the dis- 
sipative LL equation [lOj . The systems below the thresh- 
old correspond to a plane-rotator model studied in sec- 
tion 4, while the systems above the threshold are qual- 
itatively described with the present model. The men- 
tioned simulations focus on the collisions of similar-type 
(Neel or Bloch) DWs neglecting the anisotropy. They 
have predicted mutual annihilation or reflection of the 
walls depending on (parallel or anti-parallel) spin align- 
ments in both the DW centers. This result is partially 
supported by a perturbation analysis of Bloch-wall in- 
teractions within the XY model which has shown such 
DWs to repel or attract each other depending on their 
chiralities [2{| ■ The method of present study cannot be 
applied to these collisions since one is unable to determine 
neither the double-Neel nor double-Bloch wall analytical 
solutions to the dynamical equations. However, except in 
the case of a periodically distributed DWs, multi-Bloch 
or multi-Neel structures are unstable because of unbal- 
anced DW interactions, thus, they seem to be less suit- 
able for the information-storing purposes than the Neel- 
Bloch DW structures. To the best of the author knowl- 
edge, the field-induced collision of the Neel DW with the 
Bloch DW has not been simulated. 



IV. PLANE-ROTATOR MODEL 

In order to describe the DW dynamics below the 
Walker breakdown, we consider a system of plane rota- 
tors. Let us reduce the primary (LLG) dynamical system 
to its single component. Saturating the magnetization 
dynamics to the easy plain (m y — 0), we neglect the spin 
rotation about the x-axis and z-axis since the relevant 
torque components are equal to zero. For H = (H x , 0, 0), 
inserting 



m x = M 



.1 



1 + af 



m v = 0, 



= 2M- 



«i 



1 



(21) 



(where a\ takes real values) into the y-component of |T]), 
one arrives at a nonlinear diffusion equation 



dai 



dai T d 2 ai 



jH x ai - 8f3—^ + J 



dx 
2 



dx 2 



(1 



-2Ja x 



1-af) =0.(22) 



We use another ansatz describing the dynamics con- 
strained to the rry-plain (a hard plain) 



= M- 



1 



1 



= 2M- 



a 2 



1 



2 ■ 



0. (23) 



Then, we insert it into the z-component of ([T]) and we 
arrive at 



-a- 



da 2 
dt 

-2Ja 2 



^H x a 2 - 6(3^ + J 



da 2 Y 
dx J 



dx 
(0i 



r d 2 a 2 
dx 2 

■02)02(1 



(1 



0(24) 



which differs from (|22p by a constant at the anisotropy 
term. With relevance to the case 8 — 0, one finds the 
two-domain solution 



Ol(2) 



01+02 



J 



(25) 



which correspond to the Bloch DW (to the Neel DW). 
When H x ^ 0, the DW propagates with the velocity 



Cl(2) 



l\H x 



(2)1 



(26) 



The applicability of the plane-rotator model is limited 
by the Walker-breakdown condition. The magnetic field 
\H X \ cannot exceed a critical value Hyy (see Fig. 2a) 
which corresponds to the spin deviation from the basic 
magnetization plane (a canting) at the center of the DW 
about a limit angle equal or smaller than 7r/4. We esti- 
mate an upper limit of the Walker critical field consid- 
ering the x-component of the LLG equation at the DW 
center 



dm x 



dt 



El 

M 



m y m z — S 



dm x 
dx 



(27) 
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where xoi(2) = — l°g( u; )/fci(2)- Let if denotes the an- 
gle of the spin deviation (a canting) at DW center. In- 
serting ([2T|) and transforming m y m z — > m 2 z sin(<p) cos(< / s) 



in (|2"T)) or inserting 



and transforming m z m y 



dx 



(28) 



sin(i^) cos(c^) in ([27)) . one arrives at 

d°i(2) 0a . ,„ v 
—gf- = y sin(2(^)a 1(2) 

and finally, assuming \<p\ < 7r/4, at 

< H w < m&xH w = + — \k 2 \. (29) 
2 7 7 

This expression corresponds to the one given in [3(| H3, 
l38j . However, we notice that, for typical nanowires whose 
width-to-thickness ratio is bigger than 20 (double-atomic 
or triple-atomic layers of a submicrometer width), mea- 
suring Hw, one has estimated the canting angle to take 
a value of a few degrees at most, [H], f30L |39|. 

For H = 0, S ^ 0, the two-domain solution to (|2"2" ]) -(|2"4" |) 
takes the form 

a 1{2) = lue *im(»-W«). 




(30) 

It is seen that only the non-adiabatic part of the spin- 
transfer torque contributes to (|2"2")) . (|24p since the current- 
dependent term is proportional to f3. From (l28l) . the 
current-induced Walker breakdown corresponds to the 
critical current intensity 



6 w <p 2 /[2\h\(l + p/a)} 



(31) 



if H x = 0, (see Fig. 2b). It has been observed that 
Hw , Sw decrease with decreasing the nanowire width-to- 
thickness ratio, [39|,[4(|, because this ratio determines the 
strength of the easy-plain anisotropy while Hw, §w 
with /?2 — >• 0, [4l[. Notice that analytical calculations us- 
ing 2D XY model, experimental observations, and simu- 
lations of the spin ordering in nanostripes show this or- 
dering to vary along the cross-section width of the nanos- 
tripe in the DW area, thus, revealing a complex topo- 
logical structure [T3, EH- Therefore, our plane-rotator 
description is valid only for a qualitative analysis of the 
DW dynamics in the nanostripes. 

Neither finding nonstationary double-Bloch nor 
double-Neel solutions in the form of the Hirota expan- 
sion (including its second order) does not manage. In 
particular, inserting 



Ol(2) 



Wie k 1(2) x _|_ W2e k 1(2) 3 



-■yH x t/c 



Vl(2) 



\x+k 



(32) 



into (j22]) -(|24j ) . for fci 

±yj (Pi + P2)/ J, leads to the divergence of uj^) as it fol- 
lows from the approach of section 3. In order to describe 
the collision of Bloch and Neel walls below the Walker 
breakdown, I propose to apply an effective ID model as- 
suming the magnetization precession to be overdamped, 





FIG. 2. a) A scheme of the longitudinal-field dependence 
of the DW velocity for: a wire with single-axis anisotropy 
(solid line), a wire with double-axis anisotropy (dashed line), 
b) A scheme of the current-intensity dependence of the DW 
velocity for: a wire with single-axis anisotropy (solid line), a 
wire with double-axis anisotropy; ft > a (dashed line), ft < a 
(dotted line), f3 — (dash-dotted line). 



thus, taking the lhs of ([1} to be equal to zero. Inclusion 
of the constraint |m| = M leads to the modified (by ne- 
glecting the first terms on the lhs) system (H])-©. Solving 
it, we predict the field- induced DW collision below the 
Walker breakdown to result in their reflection similar to 
the one described in section 3. The reflection is accom- 
panied by the change of the Bloch wall into the Neel wall 
and vice versa. Let us emphasize that there is no spon- 
taneous DW motion below the Walker breakdown when 
neglect magnetostatic effects [43l |. 

The technological challenge of increasing the DW 
speed is especially important below the Walker break- 
down, where the driving-field is relatively weak. Refer- 
ring to this purpose, we mention an attempt utilizing an 
increase of the nanostripe-edge roughness, thus, an in- 
crease of the damping constant a, [44|. This approach 
fails since, according to simulations of |39j |. the maximum 
of the field-induced DW velocity is insensitive to a be- 
low the breakdown. It is because Cx^) °t &~ 1 \H X \ while 
\H X \ < Hw (x a. On the other hand, since /? grows with 
a, (the non-adiabatic part of the spin-transfer torque is of 
a dissipative origin), the velocity of the current- induced 
DW motion 



c = 5/3/a 



(33) 



can be insensitive to the increase of the nanostripe-edge 
roughness as well. One has attributed some reported 
DW- velocity increase due to the nanostripe-edge rough- 
ness to a decrease of its effective cross-section width. An- 
other attempt utilized an increase of Hw due to the in- 
crease of the anisotropy constant /?2- It has been done via 
the nanowire deposition on a specific crystalline substrate 
|45j. However, the most efficient method of influencing 
the maximum DW velocity below the Walker breakdown 
is the application of the transverse magnetic field, which 
is considered in the next section 46] . 
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V. DOMAIN WALL IN PERPENDICULAR TO 
EASY AXIS FIELD 

Let us define H± = H y ± iH z . For H x = 0, we search 
for a two-domain solution to ((T|) using a different 'multi- 
linearization' than used in the previous sections 

h [-iDt + JD 2 X + S(p - i)D x + aD t ] /* • g x 



The static solutions correspond to 



+ 2 f ± /r(/r/ 2 +g 2 *5i) 



Pi + y) hftsi 

- y/r 2 S 2 * = 0, 
g* 2 [-iD t - JD 2 X + 5(p - i)D x + aD t ] ft ' 9i 

(/r/ 2 + 9* 2 9i) + (pi + y) 9* 2 9ift 

+^9lh = 0, 
h9iDln ■ fl - f x 9lDl gi ■ g x = 0(34) 

Let us focus our attention on the case a,/3,8 = for 
simplity. Then one has A = / 2 = /, ffi = 92 = 9, 
while in the general case the relations ((SJ, (JTJ) apply. We 
analyze the two cases of the external-field direction; the 
one parallel to the easy plane H + = \H Z , and the one 
perpendicular to the easy plane H+ = H y . 
In the case of H+ = \H Z , we apply the ansatz 



fi = /a = 9i 
.9i = -9*2 =i(«H 



Sl e klX - llt , 

kix — lit\ 



qie 



(35) 



where k\, qi, s± denote real constants, (the parameter 
l\ can take complex values when a / 0). The solution 
in the form (|35j) describes the wall between two domains 
whose spins are deviated from the easy axis onto the 
external-field direction about an angle which grows with 
\H + \. Inserting this ansatz into (|34|) . one finds 



81 



Pi - VP! ~ l 2 Hl 



(36) 



Considering the solutions which are static in the absence 
of the electric current, li = for 6 = 0, one arrives at 



I PI -l 2 HI 



Pi J 



(37) 



In the case of H + = H y , the ansatz relating to the 
deviation of the domain magnetization from the easy axis 
onto the external-field direction takes the form 



s 2 e 



kiyX — l^t 



k 2 x—l 2 t 



fl = fl = 92 
91 = 9*2 = s 2 + 92 C" 

with real fe, 92, S2- From (|34l) . we find 



5'2 



Pl+P 2 - v /( / 3 1+/ 8 2 )2_ 7 2 ff 2 



92- 



(38) 



(39) 



l (Pi + P 2 ) 2 -l 2 H^ 

(Pl + p2)J 



(40) 



The transverse external field does not drive the DW 
motion even in the presence of the magnetic dissipation 
(a / 0). When the current through the wire and the dis- 
sipation are applied, under the transverse magnetic field, 
the DW moves with the velocity c given by (TT4")) , which is 
independent of the value of this field. Then the solution 
to (f34|) satisfies the bilinearized LLG system (Eqs. (3) 
with additional _ff+-dependent terms) at the time points 
of the discrete set t = wn/lmli, where n = 0, ±1, ±2, . . ., 
since f\ = f 2 = f, g\ = g2 = g at these points. In- 
cluding an additional to H + longitudinal component of 
the magnetic field H x drives the DW motion. For the 
realistic case Hw ~ \H X \ < \H + \ \Pi/j\, neglect- 
ing small contributions to the H^-dependent part of the 
torque, one finds the velocity of such a DW propagation 
(fT3")) or (|26[) above and below the Walker breakdown, re- 
spectively, with |fci(2)| given by (I3T1) . (|4"0|) . This velocity 
nonlinear ly increases with | H + 1, 47| . Searching for CU2), 
in the case \H X \ < Hw, additionally, I have taken the lhs 
of ((T|) equal to zero as discussed in section 4. The ma- 
nipulation of Ci(2) via the application of the transverse 
magnetic field is potentially useful for speeding up the 
processing with a magnetically-encoded information. We 
also notice that the transverse-field dependence of |fc 1 ( 2 )| 
enables influencing the magnitude of the critical current 
of the Walker breakdown 5\y , following (|3"Tj) . 



VI. CONCLUSIONS 

We have analytically studied the DW dynamics in the 
presence of the external magnetic field and the electric 
current along the magnetic wire within the LLG ap- 
proach. It has demanded overcoming the difficulty aris- 
ing from breaking the time-reversal symmetry by inclu- 
sion of the magnetic dissipation. We have removed this 
asymmetry of the dynamical system by introducing ad- 
ditional (virtual) dynamical variables, which is a similar 
trick to the Lagrangian approach to the damped har- 
monic oscillator. Determining a connection of the addi- 
tional dynamical variables to the evolution of the mag- 
netization vector in specific ranges of time, we have ana- 
lyzed the dynamics of a single DW and of a pair of DWs. 

The magnetic-field-induced velocities of the DWs, the 
formulas (fTB"]) and (|2"6"|) , and the current- induced veloci- 
ties (|T4j) and (|33|) are found to correspond to the ones 
of the Walker approach above and below the breakdown, 
respectively. According to [28j . static three-domain so- 
lutions to ID LLG equation describe pairs of Neel and 
Bloch walls. For the purposes of the qualitative dynam- 
ics analysis of a number of DWs below the Walker break- 
down, especially of the Neel-Bloch pairs, we have pro- 
posed a dynamical equation which differs from the LLG 
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one by neglecting the lhs in ([lj. Below and above the 
breakdown, the neighboring Neel and Bloch walls move 
in the presence of the longitudinal external field in the 
opposite directions. Their collision results in the DW 
reflection accompanied by the reorientation of the Neel 
wall into the Bloch wall and vice versa. In other words, 
the DWs pass trough each other without changing their 
widths and velocities, however, the head-to-head DW 
structure changes into the tail-to-tail one and vice versa. 

Our method is useful for the analysis of two-domain 
systems under the transverse (with respect to the easy 
axis) external field, which enables a verification of nu- 
merical and experimental results fl^ - li^ . A reorientation 



of the magnetic domains due to the transverse field in- 
duces a widening of the DW area up to the infinity when 
approach with the field intensity to the coercivity value. 
The consequence of the transverse-field application is an 
increase of the DW mobility (the ratio Cim/\H x \) and 
an increase of the critical current (a shift of the current- 
driven Walker breakdown). 
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